Cooperative dynamics in charge-ordered state of a-(BEDT-TTF) 2 l3 
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Electric-field-dependent pulse measurements are reported in the charge- ordered state of a-(BEDT- 
TTF)2i3. At low electric fields up to about 50 V/cm only negligible deviations from Ohmic behavior 
can be identified with no threshold field. At larger electric fields and up to about 100 V/cm a 
reproducible negative differential resistance is observed with a significant change in shape of the 
measured resistivity in time. These changes critically depend whether constant voltage or constant 
current is applied to the single crystal. At high enough electric fields the resistance displays a 
dramatic drop down to metallic values and relaxes subsequently in a single-exponential manner to 
its low-field steady-state value. We argue that such an electric-field induced negative differential 
resistance and switching to transient states are fingerprints of domain-wall dynamics inherent to 
cooperative bond-charge density wave with ferroelectric-like nature. 

PACS numbers: 71.27.+a, 71.45.-d, 71.30.+h, 71.45.-d 
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I. INTRODUCTION 

Diverse charge structures and complex nonlinear dy- 
namics in anisotropic condensed matter systems with 
electron-phonon couplings and strong Coulomb interac- 
tions have been in focus of intense scientific activity in 
recent years. Conventional charge and spin density waves 
(CDW/SDW) previously found in highly anisotropic in- 
organic and organic systems with prevailing electron- 
phonon (electron-electron) coupling feature the collective 
phason mode with well established fingerprints in applied 
ac and dc electric fields^ - — Since the phason is pinned to 
random impurities or the commensurate lattice, it can be 
observed as a narrow resonance at terahertz frequencies 
and as a broad loss peak in radio-frequency range. Other 
famous fingerprints of the phason are the nonlinear con- 
ductivity above a finite threshold field and narrow-band 
noise due to sliding. On the other hand, in the presence of 
strong onsite and intersite Coulomb repulsion a wealth of 
broken symmetry insulating ground states with spatially 
inhomogeneous charge (spin) structures is found in var- 
ious transition metal oxides and quasi-two-dimensional 
(quasi-2D) organics. 5'- In particular, organics with rich 
variety of lattice structures provide a nice playground to 
study electronic phases ranging from localized Wigner- 
crystal-like charge (spin) to delocalized CDW-like (SDW) 
modulation of charge (spin) densities. However, there are 
no well-developed theoretical models about collective ex- 
citations in these phases and how they should respond to 
applied ac and dc electric fields. 

Under applied ac electric fields, anisotropic phason- 
like dispersions were found similar to those of conven- 
tional CDWs, while at high dc electric fields contra- 
dictory results were reported ranging from negligibly 
small nonlinearities with no threshold field to the well 
defined effects due to sliding similar to conventional 
CDWs. 7 ~ 9 At much higher fields the giant nonlinear 



conductivity and negative differential resistance (NDR) 
were reported^— which were very different from those 
known in conventional CDWs. Distinct mechanisms to 
explain these remarkable nonlinear effects have been sug- 
gested ranging from melting of charge order in 0-(BEDT- 
TTF) 2 CsM'(SCN) 4 (A/' = Zn, Co) compounds^ to di- 
electric breakdown and depinning and sliding of holes 



within charge stripes in La2-xSr x NiO 
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In the 



TMET-TTP organic system the electrothermal model 
originally developed to account for the hot electron dy- 
namics was used to reproduce the measured nonlinear 
conductivity which featured NDR with no switching. 14 
On the other hand, in (9-(BEDT-TTF) 2 MZn(SCN) 4 
(M = Cs, Rb) and a-(BEDT-TTF) 2 I 3 the I - V char- 
acteristics were found to follow a power-law behavior 
I oc V a in the range of extremely low currents (lower 
than 10 nA) and low temperatures (lower than 60 K) and 
were explained within the excitonic model with a long- 
range Coulomb interaction! 15 i 16 Further, an analogy to 
SDW sliding has been invoked to explain unusual non- 
linear effects in the case of antiferromagnetic phase of 
localized spins in MDT-TS organic compound. 17 Finally, 
when the constant voltage is applied to the sample, cur- 
rent oscillations with very low fundamental frequencies 
(10-100 Hz) have been observed in the NDR region of 
several quasi-2D organics and interpreted as the switch- 
ing frequency between low and high current states in 
analogy to the well-known Gunn effect observed in two- 
valley semiconductors^ The cause of all these remark- 
able nonlinear effects, although appearing as a common 
phenomenon in many strongly correlated insulating sys- 
tems, is still unclear and the proposed interpretations are 
highly controversial. 

a-(BEDT-TTF) 2 l3 represents one of the most promi- 
nent and best studied quasi-2D organic metals with 
charge order. Synchrotron x-ray diffraction, NMR Ra- 
man and infrared vibrational measurements (Ref. [l9| and 
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Figure 1: Schematic diagrams of the pulse measurement cir- 
cuits employed for the constant-current mode (left and middle 
panels) and for the constant- voltage mode (right panel). 



references therein) demonstrate that a-(BEDT-TTF)2l3 
should be considered as the model case of the 2D bond- 
charge density wave with concomitant formation of a pat- 
tern of structural twin domains, 2 ^ These results at the 
same time indicate that the description of charge order as 
a horizontal stripe of localized charges is over-simplified 
and that the modulation of bonds in two almost diagonal 
in-plane directions needs also to be taken into account. 
Recently, on the basis of the anisotropic electrodynamic 
response it was shown that collective excitations of such 
a charge order possess a complex anisotropic dispersion 
along diagonal as well as long two crystallographic in- 
plane directions: a broad screened relaxation due to long 
wavelength phason-like excitations and a temperature- 
independent relaxation with small amplitude due to the 
motion of domain- wall pairs between twinned domains in 
the charge order texture. 21 If this interpretation is valid, 
the dynamical properties of these collective excitations 
should also manifest themselves as a characteristic non- 
linear response in high dc electric fields. Early experi- 
ments by one of us first observed nonlinear effects in dc 
conductivity at low fields^ 2 - Further nonlinear conductiv- 
ity measurements along three crystallographic directions 
have shown NDR and voltage oscillations^ The authors 
point out that the current-voltage characteristics mea- 
sured in two-probe configuration by constant current and 
constant voltage method were essentially the same and 
oscillations could be detected only for sufficiently short 
single crystals. The picture of pure horizontal stripe 
charge order was invoked to interpret these results as 
the evidence for charge order depinning and sliding. 

In this communication we attempt to clarify these 
important issues concerning the nature and origin of 
nonlincaritics in the charge-ordered ground state of a- 
(BEDT-TTF) 2 l3. We approach the charge order as a 
2D bond-charge density wave developed in the pattern 
of ferroelectric-like twin domains. Our results give evi- 
dence that long-wavelength collective excitations of this 
charge order, while possessing phason-like dispersion, 
do not show the nonlinear response in high dc electric 
fields established for conventional CDWs. Rather, short- 
wavelength excitations with domain-wall-like dispersion 
give rise to a negative differential resistance at high dc 
electric fields and reversible switching to transient high- 
conducting states. 



II. EXPERIMENTAL 

The measurements were performed on six single crys- 
tals of varying lengths (0.2-2 mm) along three in-plane 
directions: in the a and b crystallographic axes and in di- 
agonal [110] direction. The sample selected for a detailed 
study for E || a had a cross-sectional area of 0.0032 mm 2 , 
and the distance between current and voltage contacts 
was 0.37 and 0.11mm, respectively. The electric trans- 
port in the Ohmic regime was measured between room 
temperature and 20 K as described in Ref. [l9|; qualita- 
tively the same behavior was found with no hysteresis in 
the phase transition temperature region. 

At high fields the current- voltage characteristics were 
obtained by pulse measurements in the charge-ordered 
state between 100 and 135 K, a temperature range where 
systematic errors were excluded due to sample resistances 
much smaller than input impedance of instrumentation. 
Mostly a four-point configuration was used in order to 
eliminate any possible contact influence. The experi- 
ments were carried out using two measurement modes 
as shown in Fig. [TJ In the constant-current mode, a cir- 
cuit with Keithley 6221 current-source was used (Fig. [TJ 
left panel) . In order to minimize Joule heating we applied 
rectangular pulses, typically 1-5 ms wide with a settling 
time of about 5 /its and the repetition time of 100-500 ms. 
The sample resistance was measured by monitoring the 
voltage drop V at the sample and the current I through 
load resistor R s by a Tektronix TDS3014C digital storage 
oscilloscope. We also measured the total source voltage 
U. The response was typically sampled 64 times and av- 
eraged. Interestingly, out of six samples the two with a 
distance between current contacts smaller than 0.4 mm 
showed real-time voltage oscillations, which is in accord 
with previous findings ! 23 ' 24 However, the current pulses 
turned out to feature oscillations as well, indicating that 
the current source did not succeed in regulating a con- 
stant current within the pulsed The spurious origin of 
oscillations was confirmed by the experiments using an- 
other circuit for the constant-current mode in which the 
current was fixed by means of a large resistor in series 
with an HP214B pulse generator (Fig. [TJ middle panel). 
These experiments as well as experiments done in the 
constant voltage mode never did exhibit any oscillations, 
in any of its various configurations. 

In the circuit for the constant-voltage mode a Tabor 
8023 signal generator (Fig. [TJ right panel) provided volt- 
age pulses (typical width of 1-5 ms, repetition time 100- 
500 ms), and the voltage drop and current through the 
sample is again measured by an oscilloscope. 

We have additionally recorded time-resolved infrared 
reflectivity spectra, E || a and b, during the above- 
described voltage pulses. These preliminary measure- 
ments with a time resolution of 25 fj,s have been per- 
formed in the frequency range 800-4000 cm -1 using a 
Bruker Vertex 80v spectrometer in step-scan mode as a 
trigger for the voltage-pulse generator and a Hyperion 
IR Microscope with a measurement spot between volt- 
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Figure 2: (a) Low- and (b, c, d) high-electric field behavior of a-(BEDT-TTF)2l3, E || a, at a representative temperature 125 K. 
(a) Conductivity increase normalized to its Ohmic value vs. electric field, (b) Current density vs. electric field with corresponding 
resistivity vs. time plots measured in the constant-current and (c) constant- voltage mode. The current density and electric field 
are evaluated as an average of 4-5 ms. (d) Conductivity vs. current density measured in the constant-voltage mode. The full 
line is a fit to the power law a(J, T) = 0"o(T) + AJ n with ao = 0.034 f2~ 1 cm~ 1 and n = 1.8 ± 0.2. Constant-current mode data 
follow the same expression with ctq = 0.04f2 -1 cm -1 and n = 2.1 ± 0.2. 



age contacts^ 



III. RESULTS AND DISCUSSION 

All measured temperatures show qualitatively the 
same transport behavior. A typical non-linear conduc- 
tivity for fields up to 50V/cm for E || a measured in 
constant-current mode is shown in Fig. [^a). The ob- 
tained behavior is in accord with that of the constant- 
voltage mode measurement. It is evident that the thresh- 
old field for sliding cannot be clearly defined^ Qualita- 
tively similar results were obtained for the conductivity 
data along the 6-axis and diagonal [110] direction as well. 
All orientations show an almost negligible nonlinearity (a 
few percent) if compared with the standard CDW nonlin- 
ear effects which amount to at least several tens of per- 
cent at twice the threshold held. 1-3 Pulse shapes show 
only a small deviation from rectangular shape. 

The sample response drastically changes in higher 
fields [Figs, [^b) and (c)]: the current density J against 
the electric field E curve deviates from the Ohmic-like 
behavior at fields higher than 50V/cm and eventually 
shows a negative resistance behavior (NDR)jS Our re- 
sults clearly demonstrate that the observed NDR which 
sets in above the threshold field E t h is accompanied with 
a significant change in shape of the measured resistivity R 
in time t which critically depends on the applied measure- 
ment mode2£. When the constant-current mode is used, 
R-t curves show only smooth changes within pulse [Fig. 
[DJb)]; conversely when the constant- voltage mode is used, 
R-t curves show a discontinuous jump to a low-resistance 
state [Fig. [He)] and an overall drop to low fields in E-J 
plots happens at significantly higher current densities. 
Delay times become longer for fields close to thresh- 
old indicating the field-dependent formation of the con- 
nected metallic-like regions. The NDR behavior becomes 
more prominent and the threshold field _E t h increases as 
the temperature lowers (not shown). Unfortunately, due 



to a rather restricted temperature range we cannot dis- 
tinguish between conventional CDWs, for which at low 
temperatures theory predicts E t h(T) oc exp(— T/T ), 29 
and dielectric breakdown associated with the charge car- 
rier avalanche process which experimentally follows the 
E th (T) oc exp(T /T) law^ In the former case T corre- 
sponds to the strength of the pinning potential, and we 
get a value of 5 K. This is close to that found for charge 
stripes (2.5-15 K), and somewhat smaller than the values 
typically found in conventional CDWs (20-100 K). On 
the other hand, our data nicely follow the phenomcno- 
logical formula for current-dependent conductivity ex- 
pected in the presence of avalanche process of current, 
i.e., charge carriers, a(J,T) — o"o(T) + AJ n , where the 
first term is Ohmic conductivity and n in the second term 
does not depend on temperature [Fig. [5] (d)]. 

We also characterize the switching from a high to low- 
resistance state as detected in the constant-voltage mea- 
surement mode, and to test for a possible memory effects 
we have applied the following procedure (Fig. [3]). First, a 
very short (1 ms) low-field pulse was applied to measure 
the sample resistance; it was followed by an adjustable 
high-field conditioning pulse (3 ms), and in the end a 
long, low-field probing pulse was applied to track time- 
dependent recovery of sample resistance (up to 100 ms). 
Our results demonstrate that the low-resistance states 
achieved by switching are transient in nature: once the 
conditioning pulse with a voltage high enough to induce 
switching into a low-resistance state is turned off, the 
memory keeps the sample metallic-like during 2 ms and 
subsequently the sample's resistance decays back in a 
single-exponential manner to its steady-state value. The 
switching to a low-resistance state and back is a highly re- 
producible phenomenon. Additionally, after the voltage 
is applied, the switching occurs at a characteristic de- 
lay time which becomes shorter cither for larger voltage 
amplitudes at a given temperature (0.8-17 ms at 125 K) 
or at higher temperatures for a given voltage amplitude 
(not shown). 
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Figure 3: Effects of two extreme conditioning pulses with 
total source voltage U on a- (BEDT-TTF) 2I3 , E || a, at a rep- 
resentative temperature 125 K measured in constant-voltage 
mode. At sufficiently high electric fields within the condi- 
tioning pulse the sample becomes metallic and in time decays 
back to the high-resistance state. 
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Figure 4: Comparison between measured resistance vs. time 
(full lines) and an estimate (dashed lines) based on Joule heat- 
ing and the dc curve, for three representative pulses with total 
source voltage U applied to <a-(BEDT-TTF)2l3 along the a- 
axis at 125 K. 



Although such a huge (by more than a factor of 100, 
see Fig. [31 left panel) and steep reduction of resistance 
under high electric field can hardly be attributed to Joule 
heating, the decay time dependency motivated us to ver- 
ify the extent of its influence. The resistance at 125 K 
before the switching was 10kf2. The switching to low- 
resistance state of about 40 was observed with delay 
time of 3 ms upon the application of 7.34 V. The total 
Joule heat generated in the sample prior to switching was 
calculated to be less than 1 /iJ. Using the measured value 
of heat capacity (575 J/molK at 125 K) 3 -^ and the sam- 
ple volume, the upper bound of temperature rise due to 
Joule heating would be less than 1 K with the assump- 
tion of homogenous temperature. Fig. 0] shows a com- 
parison between measured resistance vs. time as opposed 
to estimates based on Joule heating and dc resistivity 
curve. A similar data analysis was done for the constant- 
current mode experiments where only a smooth decrease 
of the resistance with no delay was detected (not shown). 
The results verified that even in this measurement mode 
the resistance fell faster than the Joule heating estimate, 
which implies that heating cannot be the sole cause of 
the observed effect. Preliminary time-resolved infrared 
spectroscopy experiments also indicate that the sample 
does not undergo a transition into the high-temperature 
metallic phase as it switches to a low-resistance state 
during a voltage pulse. To elaborate, an a-(BEDT- 
TTF)2l3 heated through the metal-to-insulator transi- 
tion would show characteristic changes in the E || a, b in- 
frared spectra before and after transition^ In our time- 
resolved infrared experiment these changes were notably 
absent, indeed there are only negligible differences to the 
insulating-phase spectra through the whole duration of 
a high-field switching pulse. 26 Therefore, we are confi- 
dent that the NDR with switching and the low-resistance 
state should be an intrinsic, although transient, state of 
electron system in the charge-ordered state of a-(BEDT- 
TTF) 2 I 3 . 

Nonlinear conductivity effects similar to the ones in 
a-(BEDT-TTF)2l3 have also been observed in diverse 



quasi- ID strongly correlated inorganic and organic in- 
sulating systems. Arguments were given that the ob- 
served dynamics should be collective in nature, analo- 
gous to the case of conventional CDWs. First, the de- 
crease of Eth with increasing temperature of the form 
E th {T)/E th (0) = exp(-T/T ) has been also found in 
CDWs and ascribed to the reduction of pinning by 
thermal fluctuations of the phase of the CDW order 
parameter^ Second, the magnitude of threshold fields, 
although several orders of magnitude larger than needed 
to depin a CDW (~ 0.01-1 V/cm), is at the same time 
much smaller (again by several orders of magnitude) than 
typical breakdown fields in band insulators (10 7 V/cm). 
Third, delayed switching is an effect also found in con- 
ventional CDWs^ 2 . and it was consistently explained as 
a result of coupling among CDW domains to form a co- 
herent current-carrying state. As a rule the nonlinear 
contribution to conductivity in conventional CDWs does 
not follow the J™ power law. To the best of our knowl- 
edge, NbSe3 is the only conventional CDW system which 
shows an NDR, 33 however the effect is very small and 
dwarfed by the drop in resistivity of a-(BEDT-TTF) 2 I 3 . 

In a- (BEDT-TTF) 2 13 the values of threshold fields are 
more than one order of magnitude smaller than the ones 
found in transition-metal Mott insulators and mixed- 
stack organic systemsi 30 i 34 ' 35 Further, the delay times 
(0.01-10 ms) are 10 times longer than known from stan- 
dard CDW system (0.001-0.1 ms), 32 but more than 100 
times shorter than the ones detected in transition metal 
compounds with stripe phasesJ^ Delayed switching in 
conventional CDWs was studied in relation to sign mem- 
ory charging and was monitored as transient voltage re- 
sponses to current pulses with different histories^ In 
these experiments, as in many others, periodic voltage 
oscillations were clearly resolved. Memory effects were 
attributed to the redistribution of CDW deformations 
and the appearance of voltage oscillations to the inter- 
action of the sliding CDW with impurities. Voltage- 
mode experiments on K0.3M0O3 at very low tempera- 
tures where electron screening is absent have shown be- 
havior atypical for conventional CDWs: strong nonlin- 
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earities, delayed current switching and oscillating cur- 
rent response were detected above an extremely high 
electric field of about 100V/cm.— As far as the voltage 
oscillations are concerned our results clearly show that 
the voltage oscillations in our samples, when they ap- 
peared, are only an artifact of the measurement circuit 
and are in sharp contrast to previously published two- 
contact measurements^ We are led to conclude that a 
more systematic investigation needs to be done on single 
crystals of different geometries, from different syntheses 
and using different constant-current circuits in order to 
reconcile the existing observations. 

Therefore, we consider that delayed switching detected 
in the NDR regime as a transient current response and 
the power-law current-dependent conductivity seem to be 
a somewhat different phenomenon than the one in most 
conventional CDWs. On the other hand, the observed 
behavior is rather similar to the one of strongly corre- 
lated LaSrNiO with charge stripes^ the SrCuO Mott 
insulators,— as well as in organic systems.— In the LaSr- 
NiO extremely long delay times near threshold field were 
taken as an indication of the field-dependent formation of 
conductive regions and the onset of NDR was suggested 
to be a sign of a dielectric breakdown and concomitant 
collective motion of holes within charge stripes. The di- 
electric breakdown was also suggested in an organic con- 
ductor based on BEDT-TTF due to a resistance drop by 
a factor of 10 3 ^ while in a mixed-stack charge-transfer 
system 3 ^ the NDR was attributed to either the cooper- 
ative motion of charge defects or domain walls inherent 
to the insulating states of these systems. In all these 
cases as well as in a-(BEDT-TTF)2l3 the current-power- 
law dependent conductivity was found with the exponent 
n = 1.8 implying a self- multiplication avalanche-like pro- 
cess of charge carriers which is induced by current den- 
sity. Such a process can be also seen as a cooperative per- 
colation process among randomly distributed pre-formed 
conductive domains switched on at fields above E t h- 

In case of switching from insulating to a metallic-like 
state in a- (BEDT-TTF) 2 I 3 we have to consider the fol- 
lowing possible causes. The first one could be the de- 
struction of charge order by high electric fields as de- 
scribed above for strongly-correlated systems, or as in 
PrCaMnOi 3 - Another option are electron heating effects. 
Simple estimates indicate that both of these causes can 
be ruled out. Namely, taking the upper bound for the 
inelastic electron mean free path of about 1.35 nm in 
the temperature range 100-135 the maximum en- 
ergy provided by a field of lOOV/cm is of the order of 
13.5 jLteV ~ 0.16 K. This value is still four times smaller 
than the charge order energy gap and far below the ther- 
mal energy at measured temperatures. Dielectric break- 
down can be also ruled out since no oscillatory behavior 
of current, predicted by theories, was observed in our 
samples . 18 ' 39 It is noteworthy that the avalanche-like in- 
crease of conductivity can happen at fields which are still 
lower than needed for the onset of breakdown. Here we 
also point out a recent pump-probe experiment which 



demonstrated a photoinduced transition in a-(BEDT- 
TTF)2l3 from insulating to the metallic stated The 
initial fs-scale pulse destroys the charge order; macro- 
scopic metallic domains form and eventually relax back 
to the charge-ordered state in ps-ns times. These fast 
relaxations are in contrast to the millisecond-scale high- 
conducting state found in our experiments. 

Keeping in mind the above considerations, we suggest 
the following scenario based on the 2D bond-charge den- 
sity wave developed in a-(BEDT-TTF) 2 I 3 ^ Since we 
do not find the nonlinear behavior well known in conven- 
tional CDWs at high electric fields, we conclude that the 
long- wavelength collective excitations of this charge order 
with phason-like dispersion do not respond to the applied 
dc electric field. Rather we propose that the domain- wall 
pairs formed due to the twinned nature of the charge- 
ordered phase are responsible for the observed nonlinear- 
ities. The observed millisecond scale of response is in- 
deed too slow for phasons; more likely it corresponds to 
soliton-like excitations. In the applied dc electric field the 
coupling to the dipole moments of each unit cell breaks 
the symmetry between two opposed dipole orientations. 
At high fields the rate of formation of domain wall pairs 
strongly increases and overcomes their formation due to 
thermal excitations. Subsequently, the motion of domain 
wall pairs become increasingly correlated and thus cre- 
ates growing conduction regions until percolation pro- 
motes an NDR. It is noteworthy that a qualitatively sim- 
ilar behavior is observed in the constant current mode, 
albeit without switching and thus significantly weaker. 
This result implies that the NDR itself is independent 
of measurement-mode but achieving a prominent effect 
absolutely requires measurement conditions which allow 
unlimited current self- multiplication: an uncapped, i.e., 
voltage source. 



IV. CONCLUSION 

In summary, we demonstrate the huge negative dif- 
ferential resistance and switching to transient high- 
conduction states at high dc electric fields within molec- 
ular planes in the charge-ordered state of a-(BEDT- 
TTF) 2 l3- Conversely, at low dc electric fields only neg- 
ligibly small nonlinearities without a clear-cut thresh- 
old field are observed. We argue that this behavior of 
dc conductivity is closely related to the origin of com- 
plex dielectric response at radio-frequencies. In high-field 
regime, we suggest that the effects are accounted for by 
the dynamics of short-wavelength excitations of the two- 
dimensional bond-charge density wave, i.e., domain walls 
at interfaces of two types of ferroelectric-like domains. 
On the other hand, in low-field regime, negligibly small 
effects indicate a largely frozen dynamics of collective 
long-wavelength (phason-like) excitations whose disper- 
sion at low ac fields we have detected previously as a 
broad screened relaxation mode. Further experiments 
are needed to clarify whether such a negative differential 
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resistivity phenomenon is a unique feature of the charge 
order in a-(BEDT-TTF) 2 l3, or is common to strongly 
correlated systems close to commensurability. 
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